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Abstract 

We establish a connection between the quark model and the 1/N C expansion mass 
formulas used in the description of baryon resonances. We show that a remarkable 
compatibility exists between the two methods in the light and heavy baryon sectors. 
In particular, the band number used to classify baryons in the 1/N C expansion 
is explained by the quark model and the mass formulas for both approaches are 
consistent. 

1 Introduction 

Since pioneering work [lj in the field, the standard approach for baryon spectroscopy 
is the constituent quark model. The Hamiltonian typically contains a spin independent 
part formed of the kinetic plus the confinement energies and a spin dependent part given 
by a hyperfine interaction. The quark model results are de facto model dependent; it is 
therefore very important to develop model independent methods that can help in alterna- 
tively understanding baryon spectroscopy and support (or not) quark model assumptions. 
Apart from promising lattice QCD calculations [2], large N c QCD, or alternatively the 
1/N C expansion, offers such a method. In 1974 't Hooft generalized QCD from SU(3) 
to an arbitrary number of colors SU(iV c ) [3] and suggested a perturbative expansion in 
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1/N C , applicable to all QCD regimes. Witten has then applied the approach to baryons 
[4] and this has led to a systematic and predictive 1/N C expansion method to study static 
properties of baryons. The method is based on the discovery that, in the limit N c — ► oo, 
QCD possesses an exact contracted SU(2iVy) symmetry [5] where Nf is the number of 
flavors. This symmetry is approximate for finite N c so that corrections have to be added 
in powers of 1/N C . Notice that a baryon is a bound state of N c quarks in the large N c 
formalism. 

The 1/N C expansion has successfully been applied to ground state baryons, either light 
[61 E] or heavy [SI IS]- Its applicability to excited states is a subject of current investigations. 
The classification scheme used in the 1/N C expansion for excited states is based on the 
standard SU(6) classification as in a constituent quark model. Baryons are grouped into 
excitation bands N = 0, 1, 2,. . . , each band containing at least one SU(6) multiplet, the 
band number N being the total number of excitation quanta in a harmonic oscillator 
picture. 

The purpose of the present paper is to show that there is a compatibility between the 
quark model and the 1/N C expansion methods. It is organized as follows. We first give a 
summary of the 1/N C expansion method in Sec. [2J Then we present a relativistic quark 
model in Sec. |3] and derive analytic mass formulas from its Hamiltonian in Sec. |4J The 
comparison between the quark model and the 1/N C mass formulas is discussed in Sec. [5] 
and conclusions are drawn in Sec. [6j We point out that the results summarized hereafter 
have been previously presented in Refs. [TOj [11] for the light baryons and [12] for the heavy 
baryons. This work aims at being a pedagogical overview of these last three references. 

2 Baryons in large N c QCD 
2.1 Light nonstrange quarks 

We begin with a summary of the 1/N C expansion in the case Nf = 2, but the arguments 
are similar for any Nf. The contracted SU(2iVj) symmetry is here the group SU(4) which 
has 15 generators: The spin and isospin subgroup generators Si and T a and operators 
acting on both spin and isospin degrees of freedom denoted by Gi a (i, a = 1, 2, 3). 
The SU(4) algebra is 



In the limit N c — > oo one has [Gi a , Gjb] — > which implies the existence of a contracted 
algebra. These SU(4) generators form the building blocks of the mass operator, at least 
in the ground state band (N — 0). For orbitally excited states the generators t of SO (3), 
as well as the tensor operator also appear since the symmetry under consideration 
is extended to SU(4) ® SO (3). 

In the 1/N C expansion the mass operator M has general form 



[Si,T a ] 




[Si, Sj] 





(2) 



where the coefficients q encode the QCD dynamics and have to be determined from a fit 
to the existing data, and where the operators Oj are SU(4) ® SO (3) scalars of the form 

O - 1 Q {k) . n (k) (D 

c 

Here Og is a fc-rank tensor in SO(3) and 0§p a fc-rank tensor in SU(2)-spin, but invariant 
in SU(2)-flavor. The lower index i in the left hand side represents a specific combination. 
Each n-body operator is multiplied by an explicit factor of 1/iV"™ -1 resulting from the 
power counting rules 0], where n represents the minimum of gluon exchanges to generate 
the operator. For the ground state, one has k = 0. For excited states the k = 2 tensor 
is important. In practical applications, it is customary to include terms up to 1/N C and 
drop higher order corrections of order 1/N%. 

As an example, we show the operators used in the calculation of the masses of the 
[70, 1~] multiplet up to order 1/N C included [13] (the sum over repeated indices is implicit) 

Oi = N c 1 , 2 = ^eS\ 3 = j^T a T a , o 4 = -^s's', 

5 = ^-£^ ij G ia G ja , 6 = —i i T a G ia . (4) 

Note that although 5 and Oq carry a factor of 1/N% their matrix elements are of order 
1/N C because they contain the coherent operator G m which brings an extra factor N c . 
01 = N c 1 is the trivial operator, proportional to iV c and the only one surviving when 
N c — > oo [4j. The operators O2 (spin-orbit), O5 and 0$ are relevant for orbitally excited 
states only. All the SU(4) quadratic invariants S l S' 1 , T a T a and G %a G ta should enter the 
mass formula but they are related to each other by the operator identity [7] 

{S\ S*} + {T a , T a } + 4 {G ia , G ia } = ^N C (3N C + 4), (5) 

so one can express G ia G ta in terms of S l S l and T a T a . 

Assuming an exact SU(2)-flavor symmetry, the mass formula for the ground state band 
up to order 1/N C takes the following simple form [7] 

M^c^ + c^ + ol^j, (6) 

which means that for N = only the operators Oi and O4 (spin-spin) contribute to the 
mass. 

Among the excited states, those belonging to the N = 1 band, or equivalently to 
the [70, 1~] multiplet, have been most extensively studied, either for Nf = 2 (see e.g. 
Refs. [HI [ESJ OH [T7J US]) or for N f = 3 [Tg. The N = 2 band contains the [56', 0+], 
[56, 2 + ], [70, C + ] (£ = 0, 2), and [20, 1 + ] multiplets. There are no physical resonances 
associated to [20, 1 + ]. The few studies related to the N = 2 band concern the [56', + ] 
for Nf — 2 [2D], [56,2+] for N f = 3 [21], and [70,£+] for Nf — 2 [22], later extended to 
Nf = 3 [23J. The method has also been applied [24J to highly excited non-strange and 
strange baryons belonging to [56, 4 + ], the lowest multiplet of the N = 4 band [25j . 

The group theoretical similarity of excited symmetric states and the ground state 
makes the analysis of these states simole [2T1 [241. For mixed symmetric states, the situ- 



symmetric states to that of symmetric states. This is achieved by the decoupling of the 
baryon into an excited quark and a symmetric core of N c — 1 quarks. This procedure has 
been applied to the [70, 1"] multiplet [H CDS HS1 EH M, US] and to the [70, £+] (£ = 0, 2) 
multiplets [22J, [23] . But it has recently been shown that the decoupling is not necessary 
|13j . provided one knows the matrix elements of the SU(2iVy) generators between mixed 
symmetric states. The derivation of these matrix elements is not trivial. For SU(4) they 
have been derived by Hecht and Pang [2B] in the context of nuclear physics and adapted to 
quark physics in Ref . [13] , where it has been shown that the isospin-isospin term becomes 
as dominant in A as the spin-spin term in N resonances. 

The derivation of SU(6) matrix elements between mixed symmetric states [N c — 1, 1] 
is underway [27J. 

A detailed description of the problems raised by the standard procedure [T7] of the 
separation of a system of mixed spin- flavour symmetry [iV c — 1,1] into a symmetric core 
of N c — 1 quarks and an excited quark has been given in Refs. [28| [29] . 

2.2 Inclusion of strangeness 

For light strange baryons (Nf = 3) the mass operator in the 1/N C expansion has the 
general form 

i=l i=l 

where the operators are invariants under SU(6) transformations and the operators 
Bi explicitly break SU(3)-flavor symmetry. In the case of nonstrange baryons, only the 
operators Oj contribute, see Eq. ([2]). Therefore Bi are defined such as their expectation 
values are zero for nonstrange baryons. The coefficients di are determined from the 
experimental data including strange baryons. In Eq. ([?[) the sum over % is finite and in 
practice it containes the most dominant operators. Examples of Oi and Bi can be found 
in Refs. [HJ 1231 121] • 

Assuming that each strange quark brings the same contribution AM S to the SU(3)- 
flavor breaking terms in the mass formula, we define the total contribution of strange 
quarks as [TT] 

n a AM a = ^2diBi, (8) 

i=i 

where n s = —S is the number of strange quarks in a baryon, S being its strangeness. 

2.3 Heavy quarks 

The approximate spin-flavor symmetry for large N c baryons containing light q = {u, d, s} 
and heavy Q = {c,b} quarks is SU(6)x SU(2) C x SU(2)h, i.e. there is a separate spin 
symmetry for each heavy flavor. Over a decade ago the 1/N C expansion has been gen- 
eralized to include an expansion in 1/m.g and light quark flavor symmetry breaking [S]. 
The majority of the currently available experimental data concerning heavy baryons is 
related to ground state baryons made of one heavy and two light quarks [30] . Such heavy 
baryons, denoted as qqQ baryons, have been recently reanalyzed within the combined 
1/N C and 1/tuq expansion [9], and masses in good agreement with experiment have been 
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in Refs. [31] but much work remains to be done in this field. That is why we focus here 
on the N = band for qqQ baryons only. 

Let us first consider that SU(3)-flavor symmetry is exact. In this case the mass opera- 
tor is a flavor singlet and in the combined 1 / uiq and 1/N C expansion to order 1 / nig 
it takes the following form 

M (1) = mgl + A qq + X Q + \ qqQ . (9) 
The leading order term is uiq at all orders in the 1/N C expansion. Next we have 

A qq = Co N c l + ^-J 2 gq , and X Q = N Q ^- (c' 1 + ^ j£j , (10) 

— * — * 

where J qq is identical to the total spin S qq of the light quark pair when one deals with the 
N = band. Note that A qq contains the dynamical contribution of the light quarks and 
is independent of uiq while Xq gives I/tuq corrections. The last term, X qq Q, contains the 
heavy-quark spin-symmetry violating operator which reads 



\lQ - 2 2 Jqq ■ Jq, (11) 

where Jq is identical to the spin Sq of the heavy quark. 

The unknown coefficients cq, C2, c , c 2 , and c 2 are functions of 1/N C and of a QCD 
scale parameter A. Each coefficient has an expansion in 1/N C where the leading term (in 
dimensionless units) is of order unity and does not depend on I/tuq. Thus, without loss 
of generality, by including dimensions, one can set Co = A. The quantity A, as well as the 
other coefficients, have to be fitted to the available experimental data. In agreement with 
Ref. [S], we take 

Co = A, C2 ~ A, c ~ c 2 ~ c 2 ~ A . (12) 

The inclusion of SU(3)-flavor breaking leads to an expansion of the mass operator in 
the SU (3)- violating parameter e which contains the sing let MW, an octet and a 

27-plet M( 27 \ The last term brings contributions proportional to e 2 and we neglect it. 
For we retain its dominant contribution T 8 to order N®. Then the mass formula 

becomes 

M = M (1) +eT 8 . (13) 

The flavor breaking parameter e is governed by the mass difference m s — m (where m is 
the average of the m u and masses) and therefore is e ~ 0.2-0.3. It is measured in units 
of the chiral symmetry breaking scale parameter A x ~ 1 GeV. 

3 Quark model for baryons 
3.1 Main Hamiltonian 



The quark model used here to describe baryons aims at capturing the main physical fea- 
tures of a three-auark svstem while keeping the formalism as simole as possible in order 



F-junction confining potential, one-gluon exchange potential and quark self-energy con- 
tribution added as perturbative terms. Let us now shortly describe all these ingredients. 

A baryon, seen as a bound state of three valence quarks, can be described, at the 
dominant order, by the spinless Salpeter Hamiltonian H = Y^=i a/p? + m i + Vy, where 
rrii is the bare mass of the quark i and where Vy is the confining interaction potential. 
We use the bare mass of the quarks in the relativistic kinetic energy term as suggested by 
the field correlator method [32J, but other approaches, like Coulomb gauge QCD, rather 
favor a running constituent quark mass [33]. Although very interesting conceptually, the 
influence of this choice on the mass spectra should not be so dramatic than it could have 
been expected at the first glance: First, the bare and constituent heavy quark masses 
are nearly identical. Second, the constituent light quark masses quickly decrease at large 
momentum and become similar to the bare masses; a common limit is reached for the 
excited states. The situation is thus mainly different for low-lying nnn baryons (u and d 
quarks are commonly denoted as n), where the bare mass m n can be set equal to 0, but 
where the constituent mass is about 300 MeV [33]. However, the strength of additional 
interactions like one-gluon exchange (see next section) can be tuned in both cases and 
lead to final mass spectra which are quite similar. 

Both the flux tube model [31] and lattice QCD [33] support the Y-junction picture for 
the confining potential: A flux tube starts from each quark and the three tubes meet at 
the Torricelli (or Steiner or Fermat) point of the triangle formed by the three quarks, let 
us say the ABC triangle. This point T, located at X?, minimizes the sum of the flux tube 
lengths and leads to the following confining potential Vy = a ^2 i= I, where the 

position of quark % is denoted by and where a is the energy density of the flux tubes. If 
all the angles of ABC are less than 120°, then the Toricelli point is such that the angles 
ATB, BTC, and ATC are all equal to 120°. If the angle corresponding to an apex is 
greater than 120°, the Toricelli point is precisely at this apex. 

As xt is a complicated three-body function, it is interesting to approximate the con- 
fining potential by a more tractable form. In the following, we shall use 

3 , 

i=i 

3 

Vr = ka \xi — R , 

i=l 

where R is the position of the center of mass and k is a corrective factor [36] . The accuracy 
of the replacement (fT5i) has been checked to be very satisfactory (better than 5%) in this 
last reference provided that the appropriate scaling factor is used: k Q = 0.952 for qqq 
baryons and k\ = 0.930 for qqQ baryons. For highly excited states, the contribution of 
the configurations in which the Toricelli point is located on one of the quarks becomes 
more and more important, and one could think that the center of mass approximation ( fl~5l) 
is then wrong. But in such cases the angle made by the Toricelli point and the other two 
quarks is larger than 120° and the center of mass is consequently still close to the true 
Toricelli point. The approximation (fT5l) . although being less accurate for highly excited 
states, remains however relevant. 



(14) 
(15) 



3.2 Pert ur bat ive terms 



Besides the Hamiltonian f[M|) . other contributions are necessary to reproduce the baryon 
masses. We shall add them as perturbations to the dominant Hamiltonian (TH|) . The 
most widespread correction is a Coulomb interaction term of the form 



1 <J 



J I 



arising from one-gluon exchange processes, where Qts,ij is the strong coupling constant 
between the quarks i and j. Actually, one should deal with a running form cvs(r), but 
it would considerably increase the difficulty of the computations. Typically, we need two 
values: a = &s,qq f° r a QQ P a h an d a i = a s,qQ for a QQ pair, in the spirit of what has 
been done in a previous study describing mesons in the relativistic flux tube model [37]. 
There it was found that a±/ao ~ 0.7 describes rather well the experimental data of qq 
and Qq mesons. 

Another perturbative contribution to the mass is the quark self-energy. This is due to 
the color magnetic moment of a quark propagating through the QCD vacuum. It adds a 
negative contribution to the hadron masses [38]. The quark self-energy contribution for 
a baryon is given by 

— (17) 



where //j is the kinetic energy of the quark i, that is //j = \ \fpf + rnjj, the average being 

computed with the wave function of the unperturbed spinless Salpeter Hamiltonian (fill) . 
The factors / and 5 have been computed in quenched and unquenched lattice QCD studies; 
it seems well established that 3 < / < 4 and (1.0 < 5 < 1.3) GeV [35]. The function 
77(e) is analytically known; we refer the reader to Ref. [33] for an explicit formula. It can 
accurately be fitted by 

77(e) w 1 - (3e 2 with (3 = 2.85 for < e < 0.3, 

w ^ with 7 = 0.79 for 1.0 < e < 6.0. (18) 

e 2 

Let us note that the corrections depending on the parameter 7 appear at order 1/tuq in 
the mass formula, so they are not considered in this work. 

We finally point out that the quark model we developed in this section is spin indepen- 
dent. This neglect of the fermionic nature of the quarks is the reason why such a model 
is often called "semirelativistic" : The implicit covariance is preserved, but spin effects are 
absent. Spin dependent contributions (spin-spin, spin-orbit, etc.) typically come from 
relativistic corrections to the one-gluon exchange potential. It is useful to mention that in 
our formalism such potential terms between the quarks i and j should be of the form 



Vij oc (/ij/ij) (19) 



4 Mass formulas 



4.1 The auxiliary field method 



The comparison between the quark model and large N c mass formulas would be more 
straightforward if we could obtain analytical expressions. To this aim, the auxiliary field 
method is used in order to transform the Hamiltonian (fT4l) into an analytically solvable 
one |40|. With A = ka, we obtain 



E 



pf + m? 



3 I \H 
2 



E 

3=1 



2vj 



+ — 
2 



(20) 



The auxiliary fields, denoted as /ij and Uj, are operators, and H(fii, Uj) is equivalent to H 
up to their elimination thanks to the constraints 



d^H^Vj)] = 



Pi + mf, 



X\Sj — R\ 



(21) 



(fii) is the quark kinetic energy, and (Oi) is the energy of one flux tube, the average being 
computed with the wave function of the unperturbed spinless Salpeter Hamiltonian ffT4l) . 
The equivalence relation between Hamiltonians (|T4T) and (|20~!) is H(fii, Uj) = H. 

Although the auxiliary fields are operators, the calculations are considerably simpli- 
fied if one considers them as variational parameters. They have then to be eliminated 
by a minimization of the masses, and their extremal values /z^o and Uj t o are logically 
close to (fii) and (%) respectively [40]. This technique can give approximate results very 
close to the exact ones [41J. If the auxiliary fields are assumed to be real numbers, the 
Hamiltonian (|2"U1) reduces formally to a nonrelativistic three-body harmonic oscillator, for 
which analytical solutions can be found. A first step is to replace the quark coordinates 
Xi = {xi,X2,xs} by the Jacobi coordinates x k = {R,£,,ff} defined as [42] 

R = (nixi + n 2 x 2 + ^3X3)/ fit, with ^ = ^1 + ^2 + ^3, (22) 

— # 

and £ oc x± - x 2 , ff oc (/iifi + /i 2 x 2 ) / (a*i + /i 2 ) - ^3- 

In the case of two quarks with mass m and another with mass 7713, the mass spectrum 
of the Hamiltonian (|20|) is given by (fi\ = /x 2 = v\ = z/ 2 — v by symmetry) 



M(/i, // 3 , 1/, 1/3) = + 3/2) + w„(iV„ + 3/2) + fi + v + ^lll^ 

A A I fi 3 2/i 



??? 



// 2/i 3 



^, (23) 



where 



+ 



(24) 



//iZ/ V 2 At + ^3 V V>v 

The integers N^/ v are given by 2n^/ v + where n^/ v and are the radial and orbital 

quantum numbers relative to the variable £/ff respectively. Moreover, ^£* 2 ^ an d (ff 2 ) 

are analytically known. This eventually allows to compute ((x% — X3) 2 ) and ((x 2 — X3) 2 ), 
which are needed to know the one-gluon exchange contribution. 



The four auxiliary fields appearing in the mass formula (1231) have to be eliminated by 
solving simultaneously the four constraints 



9 M M(//, /i 3 , v, i/ 3 ) = 0, d^Mi/i, n 3 , v, z/ 3 ) = 0, 

d v M(n, /i 3 , v, u 3 ) = 0, d V3 M(n, is, z/ 3 ) = 0. (25) 

This task cannot be analytically performed in general, but solutions can fortunately be 
found in the case of light and heavy baryons. 

4.2 Light baryons 

Since we do not distinguish between the u and d quarks in our quark model and commonly 
denote them as n, there are only four possible configurations: nnn, sss, nss and snn, that 
can all be described by the mass formula ( 1231) . Important simplifications occur by setting 
m n = 0, which is a good approximation of the u and d quark bare masses. However, the 
non vanishing value for m s causes Eqs. ff25l) to have no analytical solution unless a power 
expansion in m s is performed. This is justified a priori since the strange quark is still a 
light one. After such a power expansion, the final mass formula reads [TU] 



M qqq = M + n s AM 0s (n s = 0, 1, 2, 3), 
2k aa 3fa ... m 



M = 6/i ^ - AM 0s 



a/3^0 27r/i ' ' /io 



I koaao fa ( 3 (3 



2 6V3fi 2 2tt V4/ig 5 2 



k a(N + 3) 

^o = y g • (26) 

The mass formula M qqq depends only on iV = N% + N v . The contribution of terms 
proportional to N% — N v , vanishing for n s = and 3, was found to be very weak in the 
other cases by a numerical resolution of Eqs. (|25l) . 

An important feature of the above mass formula has to be stressed: It only depends 
on N the total number of excitation quanta of the system. But, this integer is precisely 
the band number introduced in large N c QCD to classify the baryon states in a harmonic 
oscillator picture. Indeed the spinless Salpeter Hamiltonian (Tl4l) has been transformed 
into a harmonic oscillator by the auxiliary field method and it is thus natural that a such 
band number appears. The great advantage of the auxiliary field method is that it allows 
to obtain analytical mass formulas for a relativistic Hamiltonian while making explicitly 
the band number used in the large N c classification scheme to appear. The origin of N is 
thus explained by the dynamics of the three-quark system and the comparison with the 
1/N C mass formulas is therefore possible. 

4.3 Heavy baryons 

A mass formula for qqQ baryons can also be found from Eq. fT23l . An expansion in m s 
is still needed to get analytical expressions, but an expansion in 1/rriQ can also be done 



since we deal with one heavy quark. One obtains [T2] 



M qqQ = m Q + M 1 +n s AM ls + AM Q (n s = 0, 1, 2), 

fa 



AM ls 



A Mr 



4/ii - 
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12/ii 



fcia „ / 2k\a 
+ 2a r 



2iV f + 3 



iV + 3 



7T//i 




2iV € + 3 



+ 2<*i 



2k\a 



+ 



fa 



N + 3 2rr V4/if 5 



2mn 



1 



2tt/x| 

4«i 2JV,, + 3 
~3 iV + 3 



6 V 2N ( + 3 



N + 3 



^ = ^ M(iV + 3y G(iV? = ^^+3 ( v/2(iV + 3) - v / 2A^+3) • 



(27) 



At the lowest order in m s and 1/rriQ, this mass formula depends only on N. However, 
when corrections are added, the mass formula is no longer symmetric in N n and N%. 
Is it still possible to find a single quantum number? The answer is yes, provided we 
make the reasonable assumption that an excited heavy baryon will mainly "choose" the 
configuration that minimizes its mass. 

The dominant correction of order 1/mq is the term that depends on the function 
G(N,N V ). The baryon mass is lowered when G(N,N V ) is minimal, that is to say for 
N v = N. The analysis of the dominant part of the Coulomb term shows that the baryon 
mass is also lowered in this case. So it is natural to assume that the favored configuration, 
minimizing the baryon energy, is N v — N and N% = 0. It is also possible to reach the 
same conclusion by checking that an excitation of type N v will keep the baryon smaller in 
average than the corresponding excitation in N%. This is favored because of the particular 
shape of the potential, having for consequence that the more the system is small, the 
more it is light. 

As for light baryons, the quark model shows that heavy baryons can be labeled by 
a single band number N in a harmonic oscillator picture. A light diquark-heavy quark 
structure is then favored since the light quark pair will tend to remain in its ground state. 
Note that the diquark picture combined with a detailed relativistic quark model of heavy 
baryons leads to mass spectra in very good agreement with the experimental data [43J. 



4.4 Regge trajectories 

The band number N emerges from the quark model as a good classification number for 
baryons. It is now interesting to focus on the behavior of the baryon masses at large 
values of N, i.e. for highly excited states. In this limit, the formula ( J26l) gives 



M: 



qqq 



12ak (N + 3) 



24 , I6fak 

—=akoao h 6 

V3 vr 



1 + 



fak (3 



7T 



5 2 



n s m 



(2f 



Our quark model thus states that light baryons should follow Regge trajectories, that 



the baryons. The Regge slope of strange and nonstrange baryons is also predicted to be 
independent of the strangeness in the 1/N C expansion method [33]. Too few experimental 
data are unfortunately available to check this result. In the heavy baryon sector, the mass 
formula (j2"7|) with N% = and N v = N becomes at the dominant order 

(M -m Q ) 2 = 8a^(N + 3). (29) 

This model predicts Regge trajectories for heavy baryons, with a slope of 8aki/ko « 7.8a 
instead of Ylak§ « 11.4 a for light baryons. 

The Regge slope for light baryons is here given by Ylak^. However, from experiment 
we know that the Regge slopes for light baryons and light mesons are approximately equal. 
For light mesons, the exact value obtained in the framework of the flux tube model is 
2iTa, a lower value than the one obtained from formula fl28|) . This is due to the auxiliary 
field method that has been shown to overestimate the masses |45j . What can be it done 
to remove this problem is to rescale a. Let us define a such that 12aA; = 27r cr; then the 
formula (I2"5j) is able to reproduce the light baryon Regge slope for a physical value a of 
the flux tube energy density. The scaling a = TTa/(6k ) will consequently be assumed in 
the rest of this paper. 



5 Large N c QCD versus Quark Model results 
5.1 Light nonstrange baryons 

The coefficients q appearing in the 1/N C mass operator can be obtained from a fit to 
experimental data. For example, the case N = is particularly simple. Equation ([6]) 
can be applied to N and A baryons. Taking N c = 3 together with = 940 MeV for 
S = 1/2, and M A = 1232 MeV for S = 3/2, we get 

c (at=o) = 2g9 MeV ^ c ov=o) = 2g2 MeV _ ^ 

Since the spin-orbit contribution vanishes for N = 0, no information can be obtained for 
c 2 . We refer the reader to Refs. [T9], [2H [221 [21] for the determination of q at N > 0. 

In the 1/N C expansion method, the dominant term c\ N c in the mass formula (j2J) 
contains the spin-independent contribution to the baryon mass, which in a quark model 
language represents the confinement and the kinetic energy. So, it is natural to identify 
this term with the mass given by the formula (|26l) . Then, for N c = 3 we have 

Km 2tt 4vr fa 

1 9 9 1 ; 9VS 3 1 ; 

Figure [1] shows a comparison between the values of c\ obtained in the 1/N C expansion 
method and those derived from Eq. (|3T|) for various values of N. From this comparison one 
can see that the results of large N c QCD are entirely compatible with the formula (I3"T|) 
provided cr = 0.163 GeV 2 , a rather low but still acceptable value according to usual 
potential models, «o = 0.4, and / = 3.6: These are very standard values. 

Equation ffl9l implies that c 2 and C4 oc /Xq 2 ■ Therefore we expect the dependence of 
N of these coefficients to be of the form 
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Figure 1: Plot of c\ (left) and AM S (right) versus the band number N. The values 
computed in the 1/N C expansion (full circles) from a fit to experimental data are compared 
with the quark model results with a = 0.163 GeV 2 , «o = 0.4, / = 3.6, and m s = 
0.240 GeV (empty circles and dotted line to guide the eyes). No data is available for 
N = 3 in large N c studies. The large N c data are nearly indistinguishable from the quark 
model prediction in the left plot. 



We see that such a behavior is consistent with the large N c results in Fig. [2j We chose 
c\ = 208 ±60 MeV so that the point with N — 1, for which the uncertainty is minimal, is 
exactly reproduced. Let us recall that the spin-orbit term is vanishing for N = 0, so no 
large N c result is available in this case. To compute the parameter c\ a fit was performed 
on all the large N c results. In this way we have obtained c° 4 = 1062 ± 198 MeV. Note that 
c° 3> c°. This shows that the spin-spin contribution is much larger than the spin-orbit 
contribution, which justifies the neglect of the spin-orbit one in quark model studies. 



5.2 Light strange baryons 

We have first to find out the values of AM S coming from the 1/N C expansion. For N = 0, 
1, and 3, they can be found in Ref. [H], and the case N = 4 is available in Ref. [24]. 
The situation is slightly more complicated in the N = 2 band due to a larger number 
of available results. We refer the reader to Ref. [11] for a detailed discussion about the 
computation of AM S in this case. 

The mass shift due to strange quarks is given in the quark model formalism by AM 0s in 
Eq. ( |26|) . A comparison of this term with its large N c counterpart is given in Fig. [TJ where 
we used the same parameters as for light nonstrange baryons. The only new parameter is 
the strange quark mass, that we set equal to 240 MeV, a higher mass than the PDG value 
but rather common in quark model studies. One can see that the quark model predictions 
are always located within the error bars of the large N c results. Except for N = 3, whose 
large N c value would actually require further investigations, the central values of AM S 
in the large N c approach are close to the quark model results and they decrease slowly 
and monotonously with increasing N. Thus, in both approaches, one predicts a mass 
correction term due to SU(3)-flavor breaking which decreases with the excitation energy 
(or N). 
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Figure 2: Values of C2 (left) and C4 (right) versus the band number N. The values 
computed in the 1/N C expansion (full circles) from a fit to experimental data are compared 
with results from formula (1521) (empty circles and dotted line to guide the eyes). No data 
is available for N = 3 in large N c studies. 



5.3 Heavy baryons 

As mentioned previously, our present study is restricted to ground state heavy baryons 
made of one heavy and two light quarks. In the 1/N C , 1/ttiq expansion the parameters 
to be fitted are A, tuq and eA x . At the dominant order, the value of A can be extracted 
from the mass combinations [8J 

1 A 3 / 2A 2 \ 

Aq = rriQ + N C A, -(Eg + 2E£) - A Q = 2— , - S Q = - ^— -J , (33) 

resulting from the mass formula (Q. Here and below the particle label represents its mass. 
A slightly more complicated mass combination, involving light baryons as well as heavy 
ones, directly leads to ttiq, that is [9] 



-(Aq + 2^q)-- 



-(2N + 3S + A + 25) - — (4A + 3S* + 25* + Q) 
8 10 



m Q . (34) 



This mass combination gives 

m c = 1315.1 ±0.2 MeV, m h = 4641.9 ± 2.1 MeV, (35a) 

while the value 

A w 324 MeV (35b) 

ensures that the mass combinations (1331) are optimally compatible with the experimental 
values for Q = c and b. A measure of the SU(3)-flavor breaking factor is given by [S] 

2Q-A Q = ^(eA x ). (36) 
The value (eA v ) = 206 MeV leads to Sn — An = 178 MeV, which is the average value of 



The new parameters appearing in the quark model are m c , m^, k\ = 0.930, and a%. 
For the other parameters we keep the values fitted in the light baryon sector. We take 
cti = 0.7ct from the quark model study of Ref. [37J. The heavy quark masses can be 
fitted to the experimental data as follows. The quark model mass formula ff27|) is spin 
independent; it should thus be suitable to reproduce the masses of heavy baryons for 
which Jy q = 0. Namely, one expects that 



M, 



nnc\]\f = o 



A r = 2286.46 ± 0. 14 MeV, M, 



nnb\]\f = Q 



A b = 5620.2 ± 1.6 MeV. (37) 



These values are reproduced by formula (|27|) with m c = 1.252 GeV and mj = 4.612 GeV. 
It is worth mentioning that we predict M nsc \ N=0 = 2433 MeV and M nsb \ N=Q = 5767 MeV 
with these parameters. These values are very close to the experimental S c and masses 
respectively. 

We can now compare the quark model and the 1/N C , 1/toq mass formulas. On the 
one hand the mass combination (1341) leads to m c = 1315 MeV and nib = 4642 MeV. On 
the other hand, the quark model mass formula f[2Tl) is compatible with the experimental 
data provided that m c = 1252 MeV and mj = 4612 MeV. Both approaches lead to quark 
masses that differ by less than 5%. Thus they agree at the dominant order, where only 
m,Q is present. 

The other parameter involved in the large N c mass formula is A. A comparison 
of the spin independent part of the mass formulas ([9]) and (1271) leads to the following 
identification for N r = 3 
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-, with n\ = a/ kiixa /Ak Q . (38) 



According to Eqs. (fT2l) and f!35bj) one has Cq = A ~ 0.324 GeV. The quark model gives 
0.333 GeV for the expression after the second equality sign in Eq. fl38|) . which means a 
very good agreement for the QCD scale A. The terms of order 1/toq lead to the identity 



2m Q AM Q \ N=Q 
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Note that to test this relation the value of mg is not needed, like for the identity (1381) . 
The large N c parameter, A = 0.324 GeV, gives for the left hand side of (|39j) c ~ A 2 = 
0.096 GeV 2 and the quark model gives for the right hand side 0.091 GeV 2 , which is 
again a good agreement. Finally, the SU(3)-flavor breaking term is proportional to eA x ~ 
m s — m ~ m s in the mass formula f|T3|) . Equations f|T3|) . fl271) . and fl36l) lead to 
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The large N c value eA x = 0.206 GeV and the quark model estimate 0.170 GeV also 
compare satisfactorily. We point out that, except for m c and m^, all the model parameters 
are determined from theoretical arguments combined with phenomenology, or are fitted on 
light baryon masses. The comparison of our results with the 1/N C expansion coefficients 
Co, c and eA x are independent of the rriQ values. So we can say that this analysis is 



An evaluation of the coefficients C2, c' 2 , and c 2 ' through a computation of the spin 
dependent effects is out of the scope of the present approach. But at the dominant order, 
one expects from Eq. ffT9l) that C2 oc [i\ 2 and c 2 ' oc \l{ . The ratio c 2 '/c2 should thus be of 
order [i\ = 356 MeV, which is roughly in agreement with Eq. ffT2]) stating that c 2 '/c2 ~ A. 
This gives an indication that the quark model and the 1/N C expansion method would 
remain compatible if the spin-dependent effects were included. 

6 Conclusions 

We have established a connection between the quark model and the 1/N C expansion both 
for light baryons and for heavy baryons containing a heavy quark. In the latter case 
the 1/N C expansion is supplemented by an I/tuq expansion due to the heavy quark. A 
clear correspondence between the various terms appearing in the 1/N C and quark model 
mass formulas is observed, and the fitted coefficients of the 1/N C mass formulas can be 
quantitatively reproduced by the quark model. 

These results bring reliable QCD-based support in favor of the constituent quark 
model assumptions and lead to a better insight into the coefficients q encoding the QCD 
dynamics in the 1/N C mass operator. In particular, the dynamical origin of the band 
number labeling the baryons in large N c QCD is explained by the quark model. 

As an outlook, we mention two important studies that we hope to make in the future. 
First, the N — 1 baryons of qqQ type are poorly known in the 1/N C , l/mq expansion. 
They should be reconsidered and compared to the quark model. Second, the ground state 
baryons made of two heavy quarks and a light quark could be studied in a combined l/iV c , 
1/rriQ expansion-quark model approach, leading to predictions for the mass spectrum of 
these baryons. 
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